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Abstract. Dynamical networks are characterized by large complex graphs 
of interactions. We suggest a procedure of simplifying the structure of such 
graphs while preserving the spectrum of their weighted adjacency matrix. As 
the process of isospectral graph reductions maintains the spectrum of the ma- 
trix up to some known set it is possible to estimate the spectrum of the original 
matrix by considering Gershgorin-type estimates associated with the reduced 
matrix. The main result of this paper is that eigenvalue estimates improve for 
all known methods as the matrix size is reduced. Moreover, our procedure of 
isospectral graph reductions is very flexible and in particular can be used to 
obtain better eigenvalue estimates of a matrix with complex valued entries to 
whatever degree is desired. 

1. Introduction 

A simple and geometrically intuitive method for estimating the eigenvalues of 
a matrix with complex valued entries was introduced by Gershgorin in |10j . This 
method gives a result equivalent to a nonsingularity result for diagonally dominant 
matrices (see theorem 1.4 [17]) which can be traced back to earlier work done by 
Levy, Desplanques, Minkowski, and Hadamard [13j|9j[14j[TT]. Gershgorin's estimate 
was later improved on by Brauer, Brualdi, and Varga [3, 4 17 . Their results were 
similar in spirit to Gershgorin's in that each assigned to every matrix A S C nxn a 
region of the complex plane containing the spectrum of this matrix. 

For Gershgorin and Brauer we will denote these regions by T(A) and JC(A) 
respectively. For the extension of Brualdi's result given by Varga we denote the 
corresponding set by B(A). We note that the Gershgorin, Brauer, and Brualdi 
regions have the property that B(A) C K,(A) C T(A) for any complex valued 
matrix A (see [T7] for details). 

In this paper, we first extend these classical results to a larger class of square 
matrices with entries in the set W consisting of complex rational functions. We 
then use this extension to improve these original estimates. The motivation for 
considering the class of matrices with entries in W arises from the following. 

Dynamical networks are typically described by large and often complex graphs 
of interactions. In the studies of such systems it has been found that an important 
characteristic of a network's structure is the spectrum of the network's adjacency 
matrix [21 [TBI UJ IE] ■ With this in mind we note that to each matrix M with entries 
in W that there is an associated graph Get That is, the matrix M — M(G) 
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is the adjacency matrix of G. We denote by G the class of graphs with adjacency 
matrices in W nx ™. 

Using the theory developed in [5] it is possible to reduce the graph G £ G to 
another smaller graph R £ G. The main result of [B] is that the spectra of G and 
R differ at most by some known finite set. We refer to this reduction process as an 
isospectral graph reduction or simply a graph reduction of G. 

In the present paper we show that by using graph reductions one can improve 
Gershgorin, Brauer, and Brualdi-type estimates of the spectra of matrices in W nxn . 
Specifically, for M(G) £ W nx ™ the regions in the complex plane for both Gersh- 
gorin and Brauer estimates of the eigenvalues of M(G) shrink as the graph G is 
reduced (see theorems |5 . 1 1 and 5.3 for exact statements). Moreover, for Brualdi-type 
estimates we give a sufficient condition under which such estimates also improve as 
the underlying graph is reduced (see theorem |5.4[ ) . Importantly, we note that as 
square matrices with complex entries belong to W" xn then this method allows for 
improved estimates of these matrices spectra as a special case. 

We also note that, for a given graph, many graph reductions are typically possi- 
ble. Hence this process is quite flexible and can be used to systematically improve 
eigenvalue estimates. Moreover, graph reductions on a typical G £ G can be used 
to estimate the spectrum of M(G) with increasing accuracy depending on the ex- 
tent to which G is reduced. In particular, if G is reduced as much as possible the 
corresponding Gershgorin region is a finite set of points in the complex plane that 
differs from the actual spectrum of M(G) by some other finite but known set of 
points. 

Moreover, the computational procedure for improving these Gershgorin, Brauer, 
and Brualdi estimates via graph reduction does not require much effort especially 
if some special structural features of the graph are known (see section 6) . In par- 
ticular, we note that an isospectral reduction of a graph has less sets from with its 
Gershgorin and Brauer regions are built than does the unreduced graph. 

This paper is organized as follows. Section 2 defines the sets W and G as well as 
the spectrum for the matrices in W nx " or equivalently the graphs G £ G. Section 
3 extends the results of Gershgorin, Brauer, Brualdi, and Varga to the class of 
matrices with entries in W. Section 4 summarizes the theory of isospectral graph 
reductions developed in [6] and uses this to improve the eigenvalue estimates of 
section 3. Section 5 contains the main results of this paper demonstrating that our 
procedure of isospectral graph reduction gives better estimates of the spectra of 
matrices than these other methods. Section 6 gives some natural applications of 
the theorems of section 5. These include estimating the spectrum of a Laplacian 
matrix of graph, estimating the spectral radius of a matrix, and determining which 
specific reductions to use for a given matrix. 



2. Preliminaries 

In this paper we consider two equivalent mathematical objects. The first is the 
set of graphs (i.e. structures or "topologies" of networks) consisting of all finite 
weighted digraphs with or without loops having no parallel edges and edge weights 
in the set W of complex rational functions described below. We denote this class of 
graphs by G where G" is the set of graphs in G having n £ N vertices. The second 
set of objects we consider are the weighted adjacency matrices associated with the 
graphs in G. 
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By way of notation we let the weighted digraph G E G be the triple (V,E,oj) 
where V and E are the finite sets denoting the vertices and edges of G respectively, 
the edges corresponding to ordered pairs (v, w) for v, w € V. Furthermore, to : E — > 
W where u(e) is the weight of the edge e for e £ E. We will use the convention 
that w(e) = if and only if e ^ E. 

For convenience, any graph that is denoted by some triple, e.g. G = (V,E,u>), 
will be assumed to be in G. Moreover, if the vertex set of the graph G = (V, E, oj) is 
labeled V — {v\, . . . , v n } then we denote the edge (uj, Vj) by ey. For the remainder 
of this paper if G = (V, E, u>) is a graph in G n then we will assume that its vertex 
set has some labeling V = {vi, ... ,v n }. 

In order to describe the set of weights W let C[A] denote the set of polynomials 
in the single complex variable A with complex coefficients. We define the set W to 
be the set of rational functions of the form p/q where p,q £ C[A] such that p and q 
have no common factors and q is nonzero. 

The set W is then a field under addition and multiplication with the convention 
that common factors are removed when two elements are combined. That is, if 
p/q,r/s <= W then p/q + r/s = (ps + rq)/(qs) where the common factors of (ps + rq) 
and (qs) are removed. Similarly, in the product pr/qs of p/q and r/s the common 
factors of (pr) and (qs) are removed. We let C[A]" X ™ and W nx ™ denote the set of 
n x n matrices with entries in C[A] and W respectively. 

In order to stress the generality of considering the set G we note that graphs, 
which arc cither undirected, unweighted or have parallel edges, can be considered 
to be graphs in G. This is done by making an undirected graph G into a directed 
graph by orienting each of its edges in both directions. Similarly, if G is unweighted 
then it can be made weighted by giving each edge unit weight. Also multiple edges 
between two vertices of G may be considered a single edge by adding the weights of 
the multiple edges and setting this to be the weight of this single equivalent edge. 

To introduce the spectrum associated to a graph G £ G we will use the following 
notation. If G = (V, E, ui) then the matrix M(G) = M(G, A) defined entrywise by 

is the weighted adjacency matrix of G. 

We let the spectrum of a matrix A = A(X) £ W™ x ™ be the solutions including 
multiplicities of the equation 

(1) det(A(X) — XI) = 

and for the graph G we let cr(G) denote the spectrum of M(G). The spectrum of a 
matrix with entries in W is therefore a generalization of the spectrum of a matrix 
with complex entries. 

Moreover, the spectrum is a list of numbers. That is, 

a{G) = { (a h m) : 1 < i < p, a, G C, n, e N} 

where is the multiplicity of the solutions at to equation 0, p the number of 
distinct solutions, and (<7i,7ij) the elements in the list. In what follows we may 
write a list as a set with multiplicities if this is more convenient. 

As we are mainly concerned with the properties of the adjacency matrix of 
graphs in G we note, as we have previously suggested, that there is a one-to-one 
correspondence between the graphs in G™ and the matrices W nx ™. Therefore, we 
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may talk of a graph G € G™ associated with a square matrix M = M(G) in W nx ™ 
and vice- versa without ambiguity. 

3. Spectra Estimation of Graphs in G. 

In this section we extend the classical results of Gershgorin, Brauer, and Brualdi 
(see for instance |17j ) on the spectra of matrices with complex entries to matrices 
in W raxn as well as some related results. To do so we will first define the notion of 
a polynomial extension of a graph G G G. 

Definition 3.1. If G G G™ and M(G)ij = Pij/qij where Pij,qij G C[A],(?y ^ let 
Li(G, ^)—YYj=i lij for 1 < i < n. We call the graph G with adjacency matrix 

M(G,A)^{^ G ' A)M(G ' A) - ^ j , Ki tj < n 

n [L i {G,X)M{G,X)ij ~ Li(G,X)X + X i = j ~ ~ 

the polynomial extension of G. 

To justify this name note that each M(G, A)y is an element of C[A] or M(G, A) has 
complex polynomial entries. Moreover, we have the following result. 

Lemma 3.2. If G G G then a(G) C a(G). 

Proof. For G G G n note that the matrix M(G, A) — XI is given by 

mC^-XI^^f^f^ , , ^ j for !<*<„. 

3 [L i {G,X)M(G,X) ij -L i (G,X)X i = j ~ ~ 

The matrix M(G, A) — XI is then the matrix M(G, A) — A/ whose ith row has been 
multiplied by Li(X). It follows that 

n 

det (M(<5, A) - XI) = ( Y[ L t (G, A)) det (M(G, A) - XI) 

i=l 

implying a(G) C cr(G). □ 

3.1. Gershgorin- Type Regions. As previously mentioned, a well known result of 
Gershgorin's originating from |10] gives a simple method whereby the eigenvalues of 
a matrix with complex valued entries can be estimated. This result is the following 
theorem which we formulate after introducing some standard notation. 
If A G C nx " let 

n 

(2) n(A)= l A *il> l<i<n 

be the row sum of A. 

Theorem 3.3. (Gershgorin) let A G C nx ". T/ien aZZ eigenvalues of A are 
contained in the set 



T(A) = {J{z£C:\z-A tt \<r l (A)}. 
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Geometrically, Gershgorin's theorem states that all eigenvalues of a matrix A 6 
C" xn are contained in the (Gershgorin) discs of the complex plane centered at An 
of radius Vi(A) for 1 < i < n. 



In order to extend theorem 3.3 to the class of matrices W" x ™ we use the following 



adaptation of the notation given in For G G G" let 

n 

n(G, A) = l M ( G ' A )wl for 1 < i < n 

be the ith row sum of M(G). 

Theorem 3.4. Lei G € G™. Then <r(G) is contained in the set 

n 

BWr(G) = (J {A € C : |A - M(G,X) U \ < n(G,X)}. 

i=l 

Proof. Since M(G,A) G C[A]" X ™ then for a G <r(G) it follows that M(G,a) is 
a complex valued matrix. As Lemma |3.2| implies that a is an eigenvalue of the 
matrix M(G, a) then by an application of Gershgorin's theorem the inequality 
| a - M(G,a)u\ < n(G,a) holds for some 1 < % < n. Hence, a G BW r (G). □ 

Because it will be useful later in comparing different regions in the complex plane 
for G G G™ we denote 

BWrfG), = {A G C : |A - M(G, A)«| < ^(G, A)} where 1 < i < n 

and call this the ith Gershgorin-type region of G. Similarly, we call the union of 



these n sets, given in theorem 3.4 as BWr(G), to be the Gershgorin-type region of 
the graph G. 

As an illustration of theorem |3 ,4| consider the following example. Let Q G G be 
the graph with adjacency matrix 



(3) M(G) = 



A+l 1 A+1 

A 2 A A 

2A+1 1 1 

A 2 A A 

1 



As det(M(£7, A) - XI) = (-A 5 + 2A 3 + 2A 2 + 3A + 2)/(A 2 ) then one can compute 
a(Q) = {— 1, — 1, i, — i, 2}. The corresponding Gershgorin-type region 23>Vr(£7) is 
shown in figure 1. 

We note here that BWr(G) is the union of the three regions BWr(G)i, BWr(G)2, 
and BWr(G)i whose boundaries are shown in black (on line in color). Additionally, 
the interior shading of these regions reflect their intersections and the eigenvalues 
of M(G) are indicated as points. In the figures that follow we will use the same 
technique to display similar regions. 

Remark 1. From the point of view of estimating the spectum of a graph G G G 
only the union BWr(G) of the regions BWr(G)i matter. This is also true of the 
other methods presented in this paper (e.g. Brauer and Brualdi-type regions). 



3.2. Brauer Cassini-Type Regions. Following Gershgorin, Brauer was able to 
give the following eigenvalue inclusion result for matrices with complex valued en- 
tries (see [TT]L 
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FIGURE 1. The graph Q (left) and BW r (G) (right) where a(G) = 
{ — 1, —1,2, —i, i} is indicated. 

Theorem 3.5. (Brauer) Let A E C nx " where n>2. Then all eigenvalues of A 
are located in the set 

(4) K(A)= |J {zEC:\z-A l% \\z-A ]3 \<n{A)r 3 {A)}. 

l<i,j<n 

The individual regions given by {z E C : \z — Au\\z — Ajj\ < ri(A)rj(A)} 
in equation Q are known as Cassini ovals and may consists of one or two distinct 
components. Moreover, there are (™) such regions for any nxn matrix with complex 
entries. As with Gershgorin's theorem we prove an extension to Brauer's theorem 
for matrices in W" xn . 

Theorem 3.6. Let G € G n where n > 2. Then o~{G) is contained in the set 
BW K (G)= (J {AeC: |A-M(G,A)«||A-M(G,A)jj| < r^G, AJr^G, A)}. 

1<2 j'<n 

Also, BWic(G) C BWr(G). 



Proof. As in the proof of theorem 3.4 if a S er(G) then a G cr(G) and the matrix 
M(G,a) E C nxn . Brauer's theorem therefore implies that 

\a-M(G,a)u\\at - M(G,a) j:j \ < ri(G,a)rj(G, a) 

for some pair of distinct integers i and j. It then follows that, a E BW/c(G) or 
a(G) C BW K (G). 

To prove the assertion that BW K (G) C 6W r (G) let 

(5) BWjc(G)ij = {A e C : |A - M(G, X) U \\X - M(G, X) jS \ < n(G, A)r i (G, A)} 

for distinct i and j. The claim then is that BW^G)^ C BWrfG), U BW r (G) j . 
To see this assume A e BWic{G)ij or 

| A - M(G, A)«||A - M(G, \) n \ < n(G, A)r,(G, A). 

If r t (G, \)r 3 (G, A) = then cither A - M(G, A)« =_0 or A - M(G, A)^ = 0. As 
A = M(G,A),-i implies A € BW r (G) t and A = M(G, A)# implies A e 6W r (G) 3 
then A e SW r (G), U BW r (G)j. 
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Figure 2. Left: The Brauer region JC{Q) for Q in figure 1. Right: 
IC{G) C T(G). 



If n(G, A)n,(G, A) > then it follows that 

\\-M(G,\) u \\(\\-M(G,\) J3 \ 



ri(G,\) 



rj (G,X) 



< 1. 



Since at least one of the two quotients on the left must be less than or equal to 1 
then A £ BWr(G)i U BWr(G)j which verifies the claim and the result follows. □ 



We call the region BWic(G) given in theorem |3.6| the Brauer-type region of the 
graph G and the region BWic(G)ij given in ^ the ijth Brauer-type region of G. 
Using theorem |3.6| on the graph Q given in figure 1 we have the Brauer-type region 
shown in the left hand side of figure 2. On the right is a comparison between 
BW K (G) and BW r (G) where the inclusion BW K (G) C BW r (G) is demonstrated. 

3.3. Brualdi-Type Regions. In this section we extend a result of Varga [17] 
which is itself an extension of a result of Brualdi [I] relating the spectrum of a 
graph with complex weights to its cycle structure. In order to state this result we 
need the following. 

A path P in the graph G = (V, E, ui) is a sequence of distinct vertices Vx,...,v m £ 
V such that e^+i £ E for 1 < i < m — 1. In the case that the vertices vi,...,v m 
are distinct, except that v\ = v m , P is a cycle. If 7 is a cycle of G we denote it by 
its ordered set of vertices. That is, if e. iii+1 £ E for 1 < i < m — 1 and e mi £ E 
then we write this cycle as the ordered set of vertices {vi, . . . ,v m } up to cyclic 
permutation. Moreover, a cycle consisting of a single vertex is a loop. 

A strong cycle of G is a cycle {v\, . . . ,v m } such that m > 2. Furthermore, if 
Vi £ V has no strong cycle passing through it then we define its associated weak cycle 
as {v.i} irregardless of whether en £ E. For G £ G we let C S (G) and C W (G) denote 
the set of strong and weak cycles of G respectively and let C(G) = C S (G) U C W (G). 

A directed graph is strongly connected if there is a path (possibly of length 
zero) from each vertex of the graph to every other vertex. The strongly connected 
components of G = (V, E) are its maximal strongly connected subgraphs. Moreover, 
its vertex set V = {vi, . . . , v n } can always be labeled in such a way that M(G) has 
the following triangular block structure 
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M(G) 



M(Si(G)) 
* M(S 2 (G)) 







* ... * M(S ro (G)) 

where §j (G) is a strongly connected component of G and * are block matrices with 
possibly nonzero entries (see [T2] or [T7] for more details). 

As the strongly connected components of a graph are unique then for G £ G™ 
we define 

h(G,X)= W ($ e (G) A) tJ \ for I <i<n 

where i £ Nt and Ne is the set of indicies indexing the vertices in Se(G). That is, 
fi(G, A) is ri(G, A) restricted to the strongly connected component containing w,-. 

If A G C" x " then we write f^G, A) = n{A) where A = M{G). This allows us to 
state the following theorem by Varga [17] which is an extension of Brualdi's original 
theorem 0]. 

Theorem 3.7. (Varga) Let A £ C" x ". Then the eigenvalues of A are contained 
in the set 

B{A)= |J {zeC: l[\z-A zl \<l[h(A)}. 

7£C(A) UjG7 «ie7 

As with the theorems of Gershgorin and Brauer this result can also be extended 
to matrices in W nx ™ as follows. 

Theorem 3.8. Let G £ G n . Then <r(G) is contained in the set 

(6) BW B (G) = |J {A G C : J] |A-M(G,A)«| < J[ n(G,X)}. 

Also, BWb(G) C BWjc(G). 

For G G G™ we call BWb(G) the Brualdi-type region of the graph G and 
BW B (G) 7 = {A G C : [] |A-M(G,A)«|< JJ f<(G,A)} 



"iG7 



"i67 



the Brualdi-type region associated with the cycle 7. 



Proof. For G G G™, note that its polynomial extension G can be considered to be a 
graph with polynomial weights in the variable A. With this in mind let G = G(A). 
Then for fixed a G C, G(a) is the graph with adjacency matrix M(G,a) £ C" xn . 
Furthermore, since G(G(A)) is the cycle set of G(A) then for any 7 = {vx, . . . , v m } 
in G(G(A)) and fixed a £ C let 7(a) be the set of vertices {vi, . . . , v m } in the graph 
G(a). 



Using this notation, if a £ c(G) then by lemma 3.2 and theorem |3 . 7| there exists 
a 7' £ G(G(a)) such that 

(7) [] \a-M(G,a) u \< J| h(G,a). 

There are then two possibilities. Either 7' £ G(G) i.e. the set of vertices "/'(a) 
is also a cycle in G in which case a £ BWb(G) by equation (JsJ) or 7' ^ G(G). 
Suppose that the latter is the case or 7' ^ G(G). 
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Note that ifj' E C a (G(a)) then as M(G, a) t] ^ implies M(G, A)„- ^ for i ^ j 
then 7' € G S (G) which is not possible. Hence, 7' <E C^G^a)) or 7' must be a loop 
of some vertex 17 where the graph induced by {vj} in G(a) is a strongly connected 
component of G(a). Therefore, equation ^ is equivalent to \a — M(G, a)jj\ < 0. 
Hence, a — M(G,a)jj. As some cycle 7 £ G(G) contains the vertex Vj then a is 
contained in the set 

{A e C : [] |A - M(G, X) u \ < [] h(G, A)} 

implying that a € BW B (G). 

To show that SWb(G) C BW K {G) let 7 S G(G). Supposing that 7 e C W (G) 
then A = {^i} for some vertex Vi of G and 

6W B (G) 7 = {A e C : |A - M(G, A)«| = 0} 

as is the vertex set of some strongly connected component of G. It follows from 
^ that BWb{G) 1 C SWyc(G)ij for any 1 < j < n where i 7^ j. In particular, note 
that if fj(G, A) = then A S SVVx;(G)jj for any 1 < j <n where i 7^ j. 

If on the other hand, 7 € G S (G) then for convenience let 7 = . . . , v p } where 
p > 1 and note that 

p p 
(8) 6W B (G) 7 = {AeC:[|l A - M (G, A )"l < n^CG, A)}. 

i=l i=l 

Assuming < fj(G, A) for all 1 < i < p then for fixed A 6 KWb(G) 7 it follows 
by raising both sides of the inequality in ([8| to the (p — l)st power that 

rq s tt / lA-M(G,A) i ,HA-M(G,A) Jj K 

As not all the terms of the product in ^ can exceed unity then for some pair of 
indices I and k where 1 < £, k < p and i ^ k it follows that 

|A - M(G, A) fcfc ||A - M(G, A)«| < f fe (G, A)f,(G, A). 

Using the fact that f,(G, A) < ri(G,X) for all 1 < i < n we conclude that A S 
SWyc(G)A;f completing the proof. □ 

The Brualdi-type region for the graph Q with adjacency matrix ^ is shown in 
figure 3 (left) where we note that BWb{Q) Q BWk.{Q) or in this particular case 
that BW B {G) = BWk.{G) (right). 

4. ISOSPECTRAL GRAPH REDUCTIONS 

Here we present a method developed in j6] which allows for the reduction of a 
graph G G G while maintaining the graph's spectrum up to some known set. To 
do this we first introduce some definitions as well as some terminology that allow 
us to be precise in the formulation of an isospectral graph reduction. We note that 
all results in this section, with the exception of proposition 1, can be found in [6 
as well as their proofs. 

In the following if S C V where V is the vertex set of a graph we let S denote 
the complement of S in V. Also if {vi, . . . , v m } is a path in G € G let the vertices 
i»2) • • • j Vm-l °f P be its interior vertices. If P = {v\, . . . , v rn } is a cycle where we 
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FIGURE 3. Left: The Brualdi-type region BWb(G) for Q in figure 
1. Right: BW B (G) = BW k (G). 



fix some Vi € P then we say P is a cycle from Vi to Vi where P\ {vi\ are its interior 
vertices. 

Recall from section 2 that if we write the graph G as some triple (V, E, u>) then 
we are assuming Get With this in mind we give the following definitions. 

Definition 4.1. For G = (V,E,cu) let i(G) be the digraph G with all loops re- 
moved. We say the nonempty vertex set S C V is a structural set of G if 5 1 induces 
no cycles in £(G) and for each V{ € 5, w(e^) ^ A. We denote by st(G) the set of 
all structural sets of G. 

Definition 4.2. For G = (V,E,v) with S = {v u ... 1 v m } G st(G) let Bij(G;S) 
be the set of paths or cycles from vi to in G having no interior vertices in S. 
Furthermore, let 

6 S (G) = |J B l3 {G-S). 

l<i ;< 7<m 

We call the set Bs(G) the set of all branches of G with respect to S. 

Definition 4.3. Let G = (V,E,lu) and f3 e B S (G) for some S € st(G). If /3 = 
Ui, . . . , v m , m > 2 we define 

m-l , , 



A - u)(eu) 



as the branch product of /3. If m = 2 we define = w(ei2). 

is that we need Vui(P) to be 



4.1 



The reason we require w(eu) 7^ A in definition 
defined in the following. 

Definition 4.4. Let G = (V,E, u>) with structural set S — {v\ ...,v, n }. Define 
the graph IZs(G) — (S, £, n) to be the graph such that € £ if Bij(G; S) ^ and 

^])= ^2 Pw{P), l<i,j<m. 

l3eB zj (G;S) 

We call IZs(G) the isospectral reduction of G over 6 1 . 
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)A 2 -1 



fc{vi,v 2 }(-A) TZ {vi} (TZ {vuV2 }(A)) 
Figure 4. 

In order to understand the extent to which the spectrum of a graph is main- 
tained under different reductions, that is reductions of the same graph over differ- 
ent structural sets, we introduce the following. If S is a structural set of the graph 
G = {V, E, to) where V = {v x , . . . , v n } let 

N~(G;S)= [j {A € C : A = w(e«) or uj(eu) is undefined}. 

Vies 

That is, Af(G; S) is the set of A € C for which there is some vertex Vi £ S where 
w(eij) = A or, as uj(eu) — pi(X)/qi(X) £ W, the values of A at which qi(X) = 0. 

Note that, uj(en) may be zero for a given graph i.e. en could be a loop with 
weight zero and therefore not appear on the graph. However, it is important to note 
that if such is the case and Vi S S where S S st(G) then by definition £ A/"(G, S). 
That is, in any graph G = (V, E, uj) each vertex is considered to have a loop of 
possibly weight zero. 

If G, H £ G and TV is a subset of C then let cr(G) \ N be the list given by 

a{G)\N= {{a u ni) £ a(G) :a t ^N}. 

Moreover, if a(G) \ N = cr(H) \ N then we say cr(G) and cr(H) differ at most by 
the set N. The main result of [6] is the following theorem. 

Theorem 4.5. Let G £ G with S £ st(G). Then a{G) and <t(K s (G)) differ at 
most by M{G\S). 

To better understand the extent to which the spectrum of a graph G and the 
spectrum of its reduction lZg(G) differ we give the following. 

Proposition 1. Let G = (V,E,u>) be in G where S £ st(G) and S ^ V. Then 
det (M(G) - XI) 



det (M(K S (G))-\I). 



That is, the characteristic polynomial associated with the reduced graph lZs(G) 
is the characteristic polynomial of G divided by the product Y\ v . e g{^{&u) — A). As 
a consequence, for a given graph G and structural set 5* £ st(G) it is possible that 
either a(G) = a(U s {G)), a{G) C a(K s (G)), or a(TZ s (G)) C a{G). To illustrate 
this we consider the following example. 

Let A £ G be given above as in figure 4. One can compute that 

det(M(^4) - XI) = 2A - A 3 implying a(A) = {0, ±v^}. 

By proposition 1 it then follows that 

det(M(K {vuV2} {A)) - XI) = (2A - A 3 )/A 

as w(e33) = 0, since 633 is a loop of weight zero. It follows that a(lZ{ VltV2 }(A)) = 
{±V2}ca(A). 
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If the graph TZs VltV2 \(A) is reduced over the vertex {vi} then proposition 1 
implies 

A 2 — 2 2A — A 3 

det (M(K {vi} (K {vuV2} (A))) - XI) = -^-^ = j^y- 
Therefore, a(TZ {vuV2} (A)) C a(K {vl} (K {vuV2} (A))) = {0,±V2}. 



Note that for both of these reductions theorem |4.5| indicates the sets by which 
the eigenvalues of these graphs can differ. However, proposition 1 can be used to 
find not only these numbers but also the multiplicity of these numbers by which 
these spectrums can differ. 

For example, by proposition 1 the graph 1Zs V2 \(A) has two less zeros in its spec- 
trum than A provided that a(A) contains at least two zeros, otherwise a(lZ{ V2 }(A)) 
contains no zeros. As <r(-4) contains only one zero this is in fact the case. 

It is also possible for the spectrum of a graph to remain unchanged under reduc- 
tion. This happens for instance in the reduction of the graph (Jo to the graph Q\ 
where M(Q ) and M(Qi) are given on page 14 of this paper. 

Now, as any reduction IZs(G) of a graph G G G is again a graph in G it is 
natural, as in the example above, to consider sequences of reductions on a graph 
as well as to what degree a graph can be reduced. 

Definition 4.6. For G = (V,E,u) suppose the sequence of sets Si,... ,S m C V 
are such that Si G st(G), fti(G) = H Sl (G) and 

S,+i e st(Ui{G)) where n Si+1 (Ki(G)) = U l+1 (G), 1 < i < m- 1. 

If this is the case then we say Si, ... , S m induces a sequence of reductions on G 
and we write IZi(G) — 1Z(G; Si, . . . , Si) for 1 < i < m — 1. Moreover, we let 

W(G; Si, . . . , Si) = N(G; Si, ... , S,_i) U AT(7^_i(G); S ( ), 2 < i < m. 

Definition 4.7. For p £ C[A] let deg(p) be the degree of p and for ui = p/q € W let 

tt(uj) = deg(p) — deg(q). Let G™ C G n be the set of graphs where for any G G G™, 
ir(M(G)ij) < for all 1 < i, j < n. Furthermore, let G^ = U„>i G". 

Remark 2. It is important to note that any graph G where M(G) € C™ x " is a 
graph in the set G™. 

Remark 3. Note that if the sets Si, ... , S m induce a sequence of reductions on 



a graph G € G then by repeated use of theorem 4.5 it follows that a(G) and 
<t(1Z(G; Si, ... , S m )) differ at most by Af(G; Si, ... , S m ). Moreover, if the graph 
G G G,r then the set J\f(G; Si, ... , S m ) is a finite set of points in the complex plane. 

The following theorem shows that sequential reductions are in a certain sense 
commutative. 

Theorem 4.8. Let G = (V,E,oj) and V be any nonempty subset of V . If G G 
Gtt then there exists a sequence of sets Si, ... , S m _i, V C V inducing a sequence 
of reductions on G. Moreover, for any such sequence T±, . . . , T„_i, V there is a 
unique graph lZy[G] = TZ(G; T%, T n _i, V) independent of the particular sets 
T\, . . . , T m _i . 



That is, the final vertex set in a sequence of reductions completely specifics the 
reduced graph irrespective of the specific sequence of reductions. 
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5. Main Results 

In this section we give the main results of this paper. Specifically, we show that a 
reduced graph has a smaller Gershgorin and Brauer-typc regions respectively than 
the associated unreduced graph. That is, the eigenvalue estimates given in section 
3.1 and 3.2 can be improved by use of graph reductions. 

Moreover, if H is a nontrivial reduction of the graph G then there are less regions 
of the form BW r (H)i and BW K (H) i;j then of the form BW r (G), and SWjcKG) ir- 
respectively. So although there is some effort involved in reducing a graph this is 
offset by the fact that there are less Gershgorin and Brauer-type regions to compute. 

For Brualdi-type regions the situation is more complicated. For certain reduc- 
tions the Brualdi-type region of a graph may decrease in size similar to Gershgorin 
and Brauer-type regions. In other cases the Brualdi-type region of a graph may do 
the opposite and increase in size when the graph is reduced. We give an example 
of both of these possibilities in section 5.3. 

5.1. Improving Gershgorin- Type Estimates. We first consider the effect of 
reducing a graph on its associated Gershgorin region. Our main result in this 
direction is the following theorem. 

Theorem 5.1. (Improved Gershgorin Regions) Let G — (V,E,oj) where V is 
any nonempty subset ofV. //GeG, then SW r (K v [G]) Q BW r (G). 



Theorem 5.1 has the following corollary. 



Corollary 1. Let G — (V,E,to) where V is a nonempty subset of V . If the graph 
lZv[G] = 1Z{G\ Si, . . . , S m -i, V) for some sequence S±, . . . , S m -i, V then the spec- 
tmma(G)CBWr(Tl v [G})UAf(G-,S 1 ,...,S m - 1 ,V). 

In order to understand in which situations SWr (72. v [G\ ) IS strictly contained in 
BWr(G) we consider the following. For GeG" let 

dBW r (G) t = {\eC:\\-AI(G) li \ = r l (G,\)} for 1 < i < n. 

We note that the (topological) boundary of the region BWr(G)i in the complex 
plane is contained in the set dBWr(G)i. 

Theorem 5.2. Let G = (V,E,uj) be in and V be a nonempty proper subset of 
V . If Vi — V \ {vi} then at most finitely many points of the set 

U (dBW T (G)i\ |J BWr(G)j) 

are contained in ByVr(JZ-v[G]) ■ 

For a nontrivial G — (V, E, ui) we note that there is typically some region 
BWr(G)i whose boundary is not contained in the union of the other jth Ger- 
shgorin regions. Moreover, as this piece of the boundary usually contains infinitely 
many points then assuming this is the case theorem |5.1| and |5.2| imply 



BW r (K V \ {Vi} (G)) CBWr(G). 

By another application of theorem |5.1| it follows that 

BWr(n v [G}) C BWr(G) 

for any nonempty subset V C V where V does not contain Wj. That is, reducing 
over such sets strictly improve the estimates given by Gershgorin-type regions. 
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Figure 5. Left: BW r (Go)- Middle: BW r (Gi). Right: BWr(G 2 ), 
where in each the spectrum a (Go) = {~ lj — 1> — hh^} is indicated. 



As an example consider the graph Go <= with adjacency matrix 



M(G ) = 



10 1 
1 1 

10 

1 
11110 



Denoting the graphs Gi = ^{vt^.v^Go] and G2 = T^{ Vl ,v 2 }[Gi] it follows that 



M(Gi) 



A+l 

A 2 
2A+1 

A 2 





A+l 

A 
1 
A 




and M(ga) 



A+l 2A+1 



A 2 

2A+: 

A 2 



A 2 
A+l 
A 2 



The Gershgorin regions of GotGi, and G2 are shown in figure 5 where by theorem 

5^BW r (G2) C BWr(ffi) C £W r (£o)- 
Moreover, as 



dBW T (Goh \ Q BW r (5o)j and d8W r (Gi) 3 \ Q BW r(Si), 



are both infinite sets then this implies in fact that 8Wr(!?2) 5 <SWr(£i) 5 
BW r (Go)- 

Lastly, note that Af(Go\ «2) ^3}) = N{Gx] {«i, W2}) = {0}. As the point 
{0} C BWr(Si), BWr(&) then both BW r (5i) and BW T {G2) contain tr(g ) by 
corollary 1. (Note that M(Gi) = M(Q) where M(G) is given by (^J.) 

In addition an important implication of theorem |5.1| is that graph reductions on 
some G G can be used to obtain estimates of cr(G) with increasing precision 
depending on how much one is willing to reduce the graph G. Bearing this in mind 
suppose V{ is a vertex of G G G x . Then the graph 1Z{ V .}[G] = ({vi},£,{£) consists 
of a single vertex vi and possibly a loop. We note that this is the furthest extent 
to which G may be reduced. Moreover, as 7r(/i(ejj)) < 0, under the assumption 
G G G n , it follows that BWr{Tl{vi} [G]) is a finite set of points in the complex 
plane. Furthermore, if lZ{ Vi y[G] = 1Z(G; ft, . . . , S m -i, {vi}) then corollary 1 implies 
a{G) C 6W r (^ K} [G]) uV(G; ft, ... , S m -i,{vi}). As A^(G; ft, ... , S m _i,{^}) 
is also a finite set of points in the complex plane we summarize this in the following 
remark. 
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FIGURE 6. Left: BW K (Go)- Middle: BWk{Gi)- Right: BW K {g 2 ), 
where in each the spectrum a (Go) = {— 1,-1, ~i,i, 2} is indicated. 



Remark 4. Theorem |5.1| implies that graph reductions on some G £ G^ can be 
used to estimate cr(G) with increasing accuracy depending on the extent to which 
G is reduced. Moreover, if G is reduced as much as possible, i.e. to a graph on a 
single vertex, the corresponding Gershgorin-type region is a finite set of points in 
the complex plane. By corollary 1 this set differs from o~(G) by some other finite 
but known set of points. 

As an example, if G3 = TZ{ Vl y [G2] then one can compute o~(Ga) = { — 1, — lj —i, i, 2} 
and JV(0 2 ,{«i}) = {0, 1.3247, -.6623 ± 0.5622i}. Corollary 1 then implies that 
<?(Go) C {-1, -i, i, 2, 0, 1.3247, -.6623 ± 0.5622i}. 

Remark 5. The computation involved in reducing a graph G from n to m vertices 
is offset by the fact that there are n—m less ith Gershgorin-type regions to calculate. 

5.2. Improving Brauer-Type Estimates. We now consider Brauer-type regions 
for which we give results similar to those given in section 5.1. 

Theorem 5.3. (Improved Brauer Regions) Let G = (V,E,uj). If G € G^ 
where V QV contains at least two vertices, then BWjc(7Zv[G]) C BWic(G). 

Theorem |5.3| has the following corollary. 

Corollary 2. Let G = (V,E,w) where V C V contains at least two vertices. If 
the graph TZy[G] = 1Z(G; Si, ... , S m -i, V) for some sequence Si, ... , S m -i, V then 
a{G) C BW K {Tl v [G]) U N(G; Si,..., S m -i,V). 

Continuing our example, the Brauer-type regions of Go,Gi, and Gi are shown 
in figure 6 where by theorem O BW C (5 2 ) C BW K (Gi) Q BW K (Go)- More- 



over, theorem |3_6j implies C 6W r (^ 2 ), BW K (Gi) Q BW r (Gi), and 

BW K {G ) CBWr(5o). 

We note that a graph G e G must have at least two vertices for BWk.{G) to 
exist. However, a graph in G with two vertices may have a Brauer-type region 
consisting of an infinite set of points (e.g. figure 6 right hand side). That is, if 
G is reduced to a graph on two vertices, i.e. as much as is possible such that the 
Brauer-type region still exists, then this Brauer-type region will may not be a finite 
set of points. This is in contrast to the situation mentioned in remark 4 where 
Gershgorin-type regions of fully reduced graphs have this property. 

Furthermore, we note that if a graph is reduced from n to m vertices then there 
are (™) — (™) ijth Brauer-type regions to calculate. Hence, the number of regions 
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Figure 7. Left: BW B (Go)- Middle: BW B (Gi). Right: BW B (G 2 ), 
where in each the spectrum a (Go) = {— 1, — 1, — i, i, 2} is indicated. 



decrease quickly as a graph is reduced. This suggests that graph reductions may 
also decrease the amount of computation involved in estimating the spectum of a 
matrix via Brauer's method. 

5.3. Brualdi-Type Estimates. Continuing on to Brualdi-type regions we note 
that in the example we have been considering it happens that we have the inclusions 
BW B (G2) C BW b (Qi) C BW B (Go) (see figure 7). However, it is not always the 
case that reducing a graph will improve its Brualdi-type region. 

For example, consider the following graph Tt £ G w given in figure 8. If TL is 
reduced over the sets S — {i>2, «3, ttjj and T = «2> 1*3} the resulting graphs 
have adjacency matrices given by 



M(K S (H)) = 



1 

A 
10 
A 





10 


1 



and M(K T (H)) = 



1 



In this example we have the strict inclusions of the associated Brualdi-type re- 
gions given by BW B {1l s (H)) C BW B (H) C BW B (K T (H)) (see figure 8). In par- 
ticular, as BW B (lZsCH)) C BW B (H) then reducing the graph H over S increases 
the size of its Brualdi-type region. So graph reductions do not always improve 
Brualdi-type estimates of a graph's spectrum. 

However, Brualdi-type estimates still do a better job than Gershgorin and Brauer- 
type regions in estimating spectra. For instance, despite the fact that BW B (TL) 
BW B (lZs(Tt)) in the example above, theorems |3.6| and 3.8 still imply that the 



region BW B {K S {H)) C BW r (H), BW K {H), fiW r (K s (H)), BW K (Ks(H)). 

In order to give a sufficient condition under which a graph's Brualdi-type region 
shrinks as the graph is reduced we consider the following. Let G = (V,E,ui) 
where V = {vi, . . . , v n } for some n > 1 and G has strongly connected components 
§ 1 (G),...,S m (G). Define 



E" 



{e £ E : e £ S 4 (G), 1 < i < m}. 



The cycle 7 £ C(G) is said to adjacent to Vi £ V if ^ 7 and there is some 
vertex Vj £ 7 such that eji £ E scc . For any Vi £ V we denote 

A(vi, G) = {7 £ C(G) : 7 is adjacent to Vi}. 

Moreover, if C{v u G) = {7 £ C[G) : v t £ 7} then let Sfa, G) C C(v u G) be the set 
of cycles containing the following elements. 
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Figure 8. Top Left: BW B (H). Top Middle: BWb(R s (K)). Top 
Right: BWb(R-t(H)) where S — {v2,v 3 ,V4} and T= {vi,v 2 ,v 3 }. 
a(H) = {0.524, 1.490, -1.007 ± 0.513i} is indicated. 



J be a cycle in C(vi,G) where Vi 



If 



For fixed i, let 7 = {v ai , . . . 
m = 1, that is 7 = {vi}, then 7 e S(v i: G) only if 7 e C(vi,G). Otherwise, 
supposing 1 < m < n then relabel the vertices of G such that v a . is Uj for 1 < j < m 
and denote this relabelled graph by G r = (V r , £V, w r ). Then 7 e S(vi,G) if and 
only if e^i ^ S r for 1 < j < m and e m fe ^ i?^ cc for m < k < n. With this in place 
we give the following theorem. 

Theorem 5.4. (Improved Brualdi Regions) Let G = (V,E,lu) where G E 
and V contains at least two vertices. If v G V such that both A(v, G) = and 
C(v, G) = S(v, G) then BW B (n v \ v (G)) C BW B (G). 

That is, if the vertex v is adjacent to no cycle in C(G) and each cycle passing 
through v is in G) then removing this vertex improves the Brualdi-type region 
of G. We note that for graph TC in figure 8 the set A(vi,H) = {v2,vs} ^ 0. Hence, 
theorem |5.4| does not apply to the reduction of TL over S. 

However, the vertex V4 has the property that A(v4,Tt) = as well as S(v4,TC) = 
C(v4,H). Therefore, reducing H over the vertex set T = {vi, V2,i>3} improves the 
Brualdi-type region of this graph which can be seen on the upper right hand side 
of figure 8. 

As an example for why the condition C(u, G) = S(v, G) is necessary in theorem 
5.4 consider the following. Let J \H-s(J) £ G be the matrices given by 





M{J) 



then M(TZ S (J)) 



where S = {^2,^3,^4}. In this case one can compute that BWb(Ks( J)) 
BWb{J)- We note that A(v\, J) = but S{v\,J) consists of the cycle {v%, v 2 , V3} 
whereas the cycle set C{v\,J) — {{v\,v-2, U3}, {i>i, V2, i>3, ^4}}- 

Aside from this, we observe that graph reductions can increase, decrease or 
maintain the number of cycles a graph has in its cycle set. For instance the graph 
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Go in our previous example has 12 cycles in its cycle set whereas Q\ has 3 and Q 2 
has 2 (see figure 7). 

On the other hand the graphs P,TZu(P) € G with adjacency matrices given by 
10 



M{P) 



1 

1 

1 

10 10 



, and M(Ku(P)) 



10 

1 i 

1 

1 1 



where U — {vi, v 2 , v 3 , V4} have the following cycle structure. The set C(P) — 
{{vx,v 2 ,v 5 }, {v 3 ,v 3 ,v 5 }} whereas C(Ku(P)) = {{vi,v 2 }, {v 3 ,V4,}, {vi, v 2 , v 3 , u 4 }}. 
That is, reducing P over U increases the number of cycles needed to compute the 
associated Brualdi-type region from 2 to 3. 

5.4. Proofs. Here we give the proofs of the theorems in section 5.1, 5.2, and 5.3. 
The proof of theorem |5.1| is the following. 



Proof. For G — (V, E, lo) where G £ G™ and n > 2 note that the vertex set 
V\{vi} S st(G). In order to simplify notation let lZy\s Vl \(G) = TZ%, Li(G, A) = Li, 
Li(TZi, A) = L\ and M(G, X)ij — tUij. It then follows that 



BW r (Ki), = {AeC: \(X - lo v 



LOilLOu 

A - luh 



)L\\< E K w 



for 2 < i < n. Supposing A £ SWr(^i)i for fixed A and 1 < i < n then from the 
inequality above 



\(X-u, u )Ll-^LL}\< Y: \^L}\+ e . n 

By the triangle inequality 



|(A-^|-|^1^|< 
A - wn 



E K^I-kii!l + E 



i = l 



i=2 



UJilLOij 

A - W11 



ill 



A - W11 



Therefore 

(10) I (A - uu)L\\ < E \^Ll\-(\^L\\-j^\^^Ll 



i=i j/j 



j =2 



In order to use the inequality given in (10 1 let Wij = — where Pij,qij G C[A] 

Qij 



and where Lij = qa for 1 < i,j < n. Then 



(11) 



Moreover, 



BW V {G) % = {A G C : |L«(g«A - p«)| < E \Pv L v\}- 

i=i,i#j 



qiiqijqijiquX - Pu) 
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from which h\ = ]^[ qnqijqij(quX ~ pn). 



J=2 

Suppose then that 



Mi 



UaU,lj -L}\>0. 



3=2 



A - wu 



By ( 10 I this then implies that 



\(\-w«)Ll\< ]T \ u V L i\ 



As this can be written 



(12) \qa 2 (qii\-pii)L n \\L ii (q H \-p ii )\<\q? 1 2 (q n \-p 11 )L n \ ^ \PijLyl 

we consider the following. 

First, if qn = then La — implying A G BWr(G) t . Similarly, if L n = 
then A G BWr(G)i by (111. This is also the case if OnA — pn = 0. Conversely, if 
\Qif 2 (Qii^-Pii)Lii \ ^ then inequalities (11) and (12) imply that A G BW r (G) l . 
Hence, A G BW r (G)i U BW r (G) t in this case. 

On the other hand if 



(13) 



Mil 



3=2 



i^y-Lj|<o 

A - wn 



we note the following. If (<ZuA — pn)£n| = then by the above it follows that 
A G BW r (G)i U T(G) t . If |e$~ 2 (g u A -pu)Ln\ ^ then w H , and are defined 
at A for 2 < j < n and (1 1 3h may be rewritten as 



\L)\ < 



A — u>i 



J =2 



since ojh ^ in this case. Moreover, as L\/(\ — wu) cannot be zero then this 
inequality can in turn be written as |A — u>u\ < Ylj=2 l^ijl i m ply m g ^ <= BWt{G)i. 
Hence, 6W r (^i) l C BW r (G) 1 U 6W r (G) l and_ therefore BW r (^i) C BW r (G). 

By sequentially removing single vertices of V from the graph G an application 
of theorem |4.8| completes the proof. □ 



We now give a proof of theorem |5.2 



Proof. Let G G G™ for n > 1. As in the proof of theorem 5.1 let 7^y\{„ 1 }(G) = 72-i, 
Li(G, A) = Li, L<(fti, A) = L\ and M(G, A)« = w«. 
Let A G C be fixed such that 



A G <9£W r (G)i \ |J BWr(G)j 

3=2 
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It then follows that 

n 

(14) KA-wnjLil^lwyLtl 



j=2 



(15) 



| (A - Wii)it| > E 1 < i < n. 

Supposing that A G 6Wr(7?i)i for some 1 < i < n then the inequality given by 
(10 1 holds. Multiplying this inequality on both sides by yields 

n n j 

(16) il^^ka-u^i- \"i3 L i\) < ki^l(-|ii%l+£l^r^ L <l)' 



J=2 



Moreover, if both sides (14) are multiplied by |Lj/(A — wn)| then we have 
(17) 



i=2 



A - uj u 



Therefore, if l-E^-Cal ^ then (15), (16 1, and (17) yield a contradiction implying 
A^Wr(Ki),. 

If | Z/^ Z/i | = then as both Lj and L\ are polynomials in the variable A there 
are only finitely many points in the complex plane where this holds. Therefore, for 
any fixed 1 < i < n the set 

n 

(dBWr(G)i \ |J BWr(G)i) f) Bw r(^i)i 

J'=2 

consists of at most finitely many points. If V = {v\, . . . ,v m } is some nonempty, 
proper subset of the vertices of G the result follows via theorem 4^ and |5.1| by 
sequentially removing single vertices from G until V has been removed. □ 



Next we give the proof of theorem 5.3 



Proof. For G — (V, E, u>) where G £ G™ and n > 3 note that the vertex set 
V \ {vx} £ st(G). For the sake of simplicity let K V \ {vi} (G) = Tlx, L k (G,X) = 
L fc , L k (Ki,\) = L\, M(G,X) k £ = uj ke as well as R k (G) = Y%=i,&k \^uL k \ and 
A — LOkk = Afefc. 

The claim is that /Qj(fti) C /C«(G) U /Cy(G) U /Cy(G) for 2 < i, j < n where 
i j. To see this let A € K,ij{TZ\) from which it follows that 



I (A, 



)L\\\{\ 



U)j\U>ij 



) L )\ 



< 



(18) 



All An 

n n 

(Ei(^+^W)(£l(«* 



= 2 



£=2 



All 



Multiplying both sides of (18 1 by |AnLiLj| and \XnLiLj \ it follows that 
\XiiXiiLiLiL] — LViiLViiLiLiL]\\\jj\iiLiLjL^—LjjiidijLiLjL^\ < 

n n 

(E I (^^u + uau>u)LiLiLl\) ( y \ ((JjAn + Ujiuu)LiLjLj 



1=2 



£=2 
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From the triangle inequality we have 

Y\_ (\X kk XxiLiL k Ll\ — \uJkiLUikLiLkLl\j < 
ke{i,j} 



II ( \^uXxxLxL k L\\+ \u) k iu)uLiL k Ll\ ~ \uj kl uj lk Ll\) . 



ke\i,j\ 1=2 1=1 

Let A £ K, U {G), Kx 3 {G) and note then that \XxxL x \\X kk L k \ > Rx(G)R k (G) for 
k — If | An | < Rx(G) then by the previous inequality 

[] (R 1 {G)R k {G)\L\\-\u Jk xuo lk LlL k Lx\) < 

k£{i,j} 

n n 

Y\ (^Rx(G) ^ \uktL k L k \ + ^ \u k xUxkL k L k Lx\ - \uj kl uJi k L{L k Lx 

k£{i,j} 1=2 1=2 

i^k 

From the fact that 

Rx{G)R k {G)\L\\ - \u k xUx k L\L k Lx\ = 

n n 

(19) Rx(G) \"kiL l k L k \ +J2\^kxu>xkLlL k Lx\ 

t=2 1=2 

e^ k 

it follows that the previous inequality does not hold. Therefore, if A ^ ICu(G), 
tCxj(G) then \XxxLx\>Rx(G). 

Proceeding as before, where we assume again that A S tCij(lZx), then note (18) 
via the triangle inequality implies 

U {\XkkLl\-\^Ll\)< 

(20) 



ke{i,j\ i=i 11 1=2 11 

e^k 

Multiplying both sides by lAnLiijl and \XuLiLj \ furthermore implies 

^AiiLillAiiLiULj 1 ! - \uixu>xiLxLiLl\^ [\X il L l \\X j .jL j \\L 1 j \ - lujxtJxjj^LiLjL] 

(|AuLi|fli(G)|Lj| - \waUiiLiLiLll + \wixLiL\\ (Rx{G) - \XxxLx\) 

A n A 
\^XnLi\Rj{G)\L)\ - \ujxuxjLiLjLj-r^-\ + ^ {ujjxUxiLiLjLl—^l - \LOjxLjXuLi\) 

11 e=2 11 

Suppose that A £ ICxi(G), K,xj(G). Then from the above |AnLi| > Ri(G) and 
moreover if it is also the case that A £ ICij(G) then from the previous inequality it 
follows that 

» A 
Ri{G)Rj[G)\L)\ < IXuL^iG^m+y^lujjxUJxiL^L]-^] - 

1=2 Al1 
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If |AjjLj| < Ri{G) then this inequality implies 

n 

(21) lu^LjWXuLiWLh < lui^LjWXuLi] V \u u ^ 



L] 

1=2 



Note that if |AnLi| = then A G JC\k{G) for k = i,j which is not possible. Hence, 
I Anii | 7^ and multiplying both sides of (pi) by |An£i| yields 



\XnLi\\L}\ <5^|wi € i 1 ||Lj| 

implying |AnLi| < i?i(G), a contradiction. 

Therefore, if A g JCh(G), JCy(G), and /C y (G) then \X H Li\ > R^G). Moreover, 
by switching the indices i and j in this argument it similarly follows that if A ^ 
JCii(G), /Cy(G), and /C„(G) then |A i:/ Lj| > R,{G). 

Let A £ KijiJZi) and suppose that A ^ /Cii(G), /Cij(G), /Cy(G). From the above 
it follows that |AfcfcI/fc| > Rk(G) for fc = Hence, by multiplying (20 1 on both 

sides by |Lj||Lj| this implies 



n( 



R k {G)\Li\-\ U ^L k Ll\)< 



Ai 



i 



ke{i,j} 
(22) 

ke{i,j} 11 fc2 11 



Note that if 



Y l \ U ^L k Ll\-\u, kl L k Ll\>0 

An 



1=2 

and assuming as before that |An£i| ^ then by multiplying both sides of this 
inequality by |AnLi| we have 

n 

\uki L k Lj\ ^ |w w Li| > IwjtiLfci^llAn-LiI for fc = z, j. 

£=2 



As this implies that Ri(G) > |AnLi|, which is not possible, then by ( 22 1 it follows 
that A € JC U (G) U fcy(G) U /C y (G) . 

This verifies the claim which completes the proof. □ 

In order to prove theorem |5.4| we first give the following lemma. 

Lemma 5.5. Let G — (V, E, to) where G € and V contains at least two vertices. 
Moreover, suppose v% € V such that both A(v%, G) = and C(v\, G) = S(vx,G). If 
72 = {v2, ■ ■ ■ , v m } G C(TZ V \ V (G)), m > 2, and there exists a 71 = {«i, w 2 , . . . , w m } S 
G(G) ttenB>VB(^ n „(G)) 72 CB>V B (G) 

Proof. Suppose first that the hypotheses of the lemma hold for some G S G™. 
We then make the observation that the edges belonging to no strongly connected 
component of G are not used to calculate to BWb(G). Furthermore, any cycle of 
G is contained in exactly one strongly connected component of this graph. This 
implies that the Brualdi-type region of the graph is the union of the Brualdi-type 
regions of its strongly connected components. Therefore, we may without loss in 
generality assume that G consists of a single strongly connected component. 
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As in the previous proofs let lZy\{ Vl }(G) = ^l; L k (G 7 A) = Lfc, L k {lZ\, A) = L\, 
M(G,X)u = Lou as well as Rk(G) = Y^e=i,i^k \ u keLk\ and A - u kk = X kk . In 
addition, suppose that both 71 = {vx,...,v m } and 6 = {vi,v m } are cycles in 
C(vi,G) for some 1 < m < n. The fact that 71 S G(vi,G) implies in particular 
that 72 = {i>2, • ■ ■ , v m } is a cycle in C{1Z\). 

From the assumption that v\ has no adjacent cycles it follows that uj m i — for 
1 < i < m since otherwise {vi, Uj+i) . . . , v m } € -4(«i, G). Also, as 71 <E G(i>i, G) = 
S(v%, G) then w mi = for m < i < n and Uji = for 1 < j < m. Therefore, 

m m— 1 

(23) 6W B (G) 71 - {A e C : J] |A«L<| < \cj ml L m \ i? 4 (G)}, 



i=l i=l 



(24) EWb{G) 5 = {X eC:|A u Li||A im| < |w TO iL m |i2i(G)}. 
In addition, the region BWb(R-i) 12 is given by the set 

m— 1 m— 1 m— 1 

(25) {A e C : |(A mm - ^f^)^| [] |A^| < £ l"^^! II ^( G »- 

11 i=2 i=2 11 i=2 

In this case L\ = for 1 < i < in since for each such i the edge en ^ E. 
Letting u>ij = Pij/qij for all 1 < i, j ' < n then 

n 

L ln = (q m i(quX -pii))™ -1 Y[qu- 

If <?nA — pu = or q m i — then |An| = 0. Similarly q a — implies that \Xa\ = 
for all 2 < i < n. Therefore, if it is the case that = then A € £W_b(G) 7i . 

Suppose then that both A € SWb(7£i) 72 and 7^ 0. By multiplying both 
sides of (251 by |AnLi| and dividing out |£^J then by the triangle inequality 



aim i) n i a «^i ^ 

i=2 
771— 1 

\u m iL m \(Ri(G) — |wi m Li 

I) LI R*(G). 



Note that if UTJ2 1 R *(G) = then A e £W B (G) 7l and if |AnLiA mm L„ 



\u m iL x u lm L m \ = or lAuiiA^Ljnl = then A e BW B (G) S - Otherwise, ( [26] ) 
implies 

/ 2? x IE 1 |Aii-^i| < |^ m i£ m |(^i(G) - \uj lm L!\) 

UT^ 1 R i(G) ~ |AnLiA mm L m | - |w m iLiwi m L m | 
With this in mind, we note that if 

|(-Ri(G) - |cji m Za|) < |w m iL m |-Ri(G) 



I XnLiX mm L m \ — \ui m iLiLUi m L m \ \XnLiX mm L m \ 

then it follows from ( [27] ) that A S 6Ws(G) 7l . On the other hand if this inequality 
does not hold then |AnLi||A mm L m | < \w m iL m \Ri(G) or A S 23Wb(G),s. Hence, if 
<5 e G(G) then 23WbOKvv,(G)) 72 C £W B (G) 7l U SWjj(G) { . 

Conversely, if 5 <£ G(G) then |w lm L m = 0|. From ([26} it follows that BW B (G) 72 = 
BWb(G) 71 . This completes the proof. □ 
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We now give a proof of theorem |5.4| 

Proof. First suppose that the graph G = (V, E, us) in G™ consists of a single strongly 
connected component that is nontrivial. Moreover, for the vertex v\ G V suppose 
both A(vx,G) = and C(v u G) = S(v u G). 

Using the notation TZi = 7Z v \{ Vl }(G) let 7 = {v 2 , . . . ,v m } be a cycle in C(lZi) 
for some 1 < m < n. If 7 G C(G) then as A(vi,G) = it follows that M(G, A)y = 
M{Tl ll \) lj for all 2 < i < m and 1 < j < n. Hence, BW B (R-ih = 6W B (G) 7 or 
SW B (^i) 7 C BWb(G). 

On the other hand if 7 ^ C(G) then for at least one 1 < i < m both of the edges 
ii ej. « € -E- Suppose then that this is the case for only one such < i < m 
and without loss in generality that this happens at i = 1. Therefore, the cycle 
6 = {vx,.. .,v m } G C(G). Moreover, as C(vi,G) = S(vi,G) and A(vi,G) = by 



assumption then lemma 5.5 implies that EWsfKi), C BWb(G). 

If there are multiple < i < m such that both en, e-x^+x G E then we define the 
following. Denoting vq = v m then let T = {i : e,--!^, € E, 1 < i < m}. We give 
the set X the ordering X = {ii, . . . , z^} such that ij < u- if and only if j < k. Then 
the sets 

7j = { v k G 7 : ij < k < ij+i} for 1 < j < I 

along with 7^ = {vk G 7 : t>fc < ?i or > i{\ partition the vertices of the cycle 7. 

Moreover, let 7^ = jj U {«i} be given the ordering {v\, Vi j , v^+i, . . . , 
for 1 < j < I and 7] = 7, U {v x } the ordering {t>i, 1^,^+1, . ..,»*,%.. • ,^-1}. 
Then each 7J G G(t>i,G) for 1 < j < £. Furthermore, let 7^ be given this same 
ordering as jj for I < j < t with the vertex v\ removed. If this is the case then 
each 7j G C(1Zi). We note that by another application of lemma 
BW B {'R- 1 ) 1] C BW b (G) for each 1 < j < t 

The claim then is that the region 



5.5 



it follows that 



(28) BWsin^ C |J SW B (Ki) 



73 • 



To see this denote Aj< = (A - w« and i?£ = V |M(7£i, A)»j|. Then 

An . r-f,. 



. 1 , , Wii^ii , 

Vi — l A _ ' 

using this notation we have 

rn m 

BWbCJZ,), = {A G C : [] 141 < II R ^ and 



i=2 i=1 



= {A G C : J] 141 < II R i} for 1 < J < ^ 

Therefore, assuming that A ^ SW B (7^i) 7j . for all 1 < j < £ implies A ^ BW B (Ki) r 
This follows from the fact that the vertices of the cycles jj for 1 < j < t partition 
the vertices of 7. Hence, (28 1 holds or BWb(T^i)j is contained in the union of 



the regions BW bO^-i)^ which are themselves contained in BWb(G). Therefore, 
#W B (fti) 7 C BW B {G) for any 7 G C(^i). □ 



To end this section we give a proof of proposition 1 . 
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Proof. Assuming the hypotheses of the proposition suppose that S = {«i} and note 
that this in particular implies w(en) 7^ A. Let di n (vi) and d out (yi) be the in and 
out degree of v\ respectively. Using the notion of a branch expansion from [6 J let X 
be a branch expansion of G with respect to S. Then there are n = di n (yi) ■ d ou t(v\) 
branches in Bs(X) of length two each with a single intermediate vertex having a 
loop of weight w(en). Moreover, by repeated use of the proof of lemma 4.4 in [6] 
we have 

(w(en) - \) n l det (M(G) - XI) = det (M(X) - Xl) 
as there are n — 1 more vertices in X than in G each adding the solutions of 
A = w(en) to cr(G). 

Contracting the n branches of length two in £>s(G) to single edges yields 

det (M(X) - XI) = (w(e u ) - A)" det (M(TZ S (G)) - XI) 

by use of the proof of lemma 4.5 in |BJ. Therefore, the result follows in this case. 
For general structural sets S € st(G) the result follows by sequentially reducing 



G over single vertices of S by use of theorem 4.8 □ 



6. Some Applications 

In this section we discuss some natural applications of using graph reductions to 
improve estimates of the spectra of certain graphs or equivalently matrices in W" xn . 
Our first application deals with estimating the spectra of the Laplacian matrix of 
a given graph. Following this we give an algorithm for estimating the spectral 
radius of a matrix using graph reductions. Last, we use the results of theorem 
\5.2\ as well as some structural knowledge of a graph to identify particularly useful 
structural sets. The motivation for finding such sets in general is that reducing over 
them allows for better eigenvalue estimates with minimal effort. This is especially 
useful for establishing eigenvalue estimates for large graphs with known structural 
properties. 



6.1. Laplacian Matrices. An important application of theorem 4.5 is that one 
may reduce not only the graph G but also the graphs associated with both the 
combinatorial Laplacian matrix and the normalized Laplacian matrix of G. Such 
matrices are typically defined for undirected graphs without loops or weights but 
this definition can be extended to graphs in G (see remark 6 below). However, 
here we give the standard definitions as these are of interest in their own right (see 

mm- 

Let G = (V, E) be an unweighted undirected graph without loops, i.e. a simple 
graph. If G has vertex set V = {v\, . . . , v n } and d(vi) is the degree of vertex Uj then 
its combinatorial Laplacian matrix Ml(G) is given by 

!d(vi) if i = j 
— 1 if i 7^ j and Vi is adjacent to Vj 
otherwise 

On the other hand the normalized Laplacian matrix Mc(G) of G is defined as 

1 if i — J an d d(vj ) ^ 

Mc(G)ij = <j ^ d( :' )d{Vj) ^ V ^ adjacent to Vj 
otherwise 
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Figure 9. Left: BW r (L(H)). Right: BWr(n s (L(H))) , where 
in each the spectrum a(L(H)) = {0, 1, 2, 4, 5} is indicated. 



The interest in the eigenvalues of Ml(G) is that <j(Ml(G)) gives structural 
information about G (see |7j). On the other hand knowing a(M^{G)) is useful in 
determining the behavior of algorithms on the graph G among other things (see 
0). 

Let L(G) be the graph with adjacency matrix Ml(G) and similarly let C(G) be 
the graph with adjacency matrix Mc(G). Since both L(G),£(G) £ either may 
be reduced over any subset of their respective vertex sets. 

For example if H £ G v is the simple graph with adjacency matrix 



M(H) 
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1 


1 
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then the graph L(H), has the structural set S 
set yields 1Zs(L(H)) where 



{vx, V2,v 3 ,V4 : }. Reducing over this 



M(K S (L(H))) = 



r A-3 


1 


1 


1 


A-4 


A-4 


A-4 


A-4 


1 


2A-7 
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-A+5 


A-4 


A-4 


A-4 
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2A-7 


-A + 5 


A-4 


A-4 
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1 


-A+5 


-A+5 


3A-11 


L A-4 


A-4 


A-4 


A-4 



Figure 9 shows the Gershgorin regions for L(H) as well as TZs(L(H)). 

Note that the adjacency matrix of H is symmetric so its eigenvalues must be 
real numbers. With this in mind the Gershgorin- type region associated with simple 
graphs and their reductions can be reduced to intervals of the real number line. 

Remark 6. It is possible to generalize Ml(G) to any G £ G if G has no loops and n 
vertices by setting M L (G) l3 = -M{G) tJ for i ^ j and M L (G) U = Y%=ij& M ( G h- 
This generalizes and is consistent with what is done for weighted digraphs in |18j 
for example. 
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6.2. Estimating the Spectral Radius of a Matrix. For G € G p i the spectral 
radius of G denoted p(G) is the maximum among the absolute values of the elements 
in a(G) i.e. 

p{G) = max |A|. 

AScr(G) 

For many graphs G G it is possible to find some structural set S G st(G) such 
that each vertex of S has no loop. If S is such a set then a(G) and a(lZs(G)) differ 
at most by J\f(G;S) — {0}. Moreover, as Af(G;S) is in some sense the error in 
estimating a(G) by a(K s (G)) then it follows that p(G) = p(K s (G)). 

For example consider the graph K shown in figure 10 and note that the ver- 
tices Vi,V2,V3 are the vertices of K without loops. As {fi, 1*2,^3} € st{K) then 
M{G; {v\, V2, V3}) ~ {0} where lZ{ VlyV2tV3 y(G) is the graph shown on the lower 
right of the figure. 

By employing the region BWr(K) it is possible to estimate p(K) < 3. However, 
via B>Vr(K{„ h „ 2i „ 3 }(if)) it follows that p(K) < 2 (see the top left and right of figure 
10). The important idea here is that each vertex without loop may be removed from 
a graph without effecting its spectral radius. 

It should be noted that for a given graph there is often no unique set of vertices 
without loops which is simultaneously a structural set. Therefore, there may be 
many ways to reduce a graph without effecting its spectral radius. 

Moreover, it is possible to continue reducing the graph even if, as in figure 10, the 
graph has been reduced to a graph in which each vertex has a loop. For example, 
the graph K = lZ{ Vl V2 V3 y(K) in this figure can be further reduced over the set 
{vi,v 3 } where one can compute N(K; {v\, V3}) = {(1 ± v / 5)/2j. 

Since (l±V5)/2 6 BW r (K) then it follows from theorem [its] that 

a{G) C BW r (K {vi , V3} (K)) U {0, (1 ± VE)/2}. 
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Figure 11. Left: The graph N. Right: BW r (JV). 



Therefore, p(K) can be estimated using BW r (K{„ I ,„ 3 }(A')) U {(1 ± VE)/2}. 

Had it been the case that some k E JV(K; {v\, v^}) was not contained in BWr(K) 
then by our previous calculations k € cr(K) and we could then ignore it. This 
process can then be continued in the same way over some further structural set of 
TZs Vl>V3 \(K) if desired. 

To summarize this method, the first step in estimating the spectral radius of a 
graph G £ G via graph reductions is to remove those vertices of G having no loop. 
This can be accomplished either by a single reduction or sequence of reductions 
and results in a graph Q where p(G) — p(G). That is, p(G) can be estimated via 
BW r (G) rather than BW r (G). 

If greater accuracy in estimating p{G) is desired it is possible to improve this 
estimate by further reducing Q over some S £ st(Q). This is done by using the set 
BWr(K s (G)) U (N{G\ S) n BW r (G)) to estimate p{G). This is possible since 

a(G) c BW r (n s (G)) U (M(G; S) n BW T {G)) c BW r (G). 

This second step can then be repeated by reducing lZs{G) over some structural set 
and so on if improved estimates are desired. 

6.3. Targeting Specific Structural Sets of Graphs (and Networks). In this 
section we consider reducing graphs over specific structural sets in order to improve 
eigenvalue estimates when some structural feature of the graph is known. To do so 
consider G = (V, E, uu) where V = {vi, . . . , v n }. 

If the sets BWr(G)i for 1 < i < n are known or can be estimated by some 
structural knowledge of G then it is possible to make decisions on which structural 
sets to reduce over based on these sets. That is, it may be possible to identify 
structural sets V C V such that 



(J (dBW r (G)i \ |J BWr(G). 



is some infinite set. If this can be done, theorem 5.2 then implies that a strictly 
better estimate of a(G) can be achieved by reducing over V. 

For example consider the graph N = (V,E,uj) in the left hand side of figure 11 
where V — {v±, ... ,v n } for some n > 5. If it is known that N is a simple graph 
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Figure 12. Left: ft v \ {t)l} (iV). Right: i3W r (^ y \ {t , l} (iV)). 



such that d(vi) = 4, = d(v 3 ) — d{vi) — d(v$) — 3 and d{v%) G {0, 1,2,3} for 

all 6 < % < n then the sets BWr(G)i are each discs of radius either 0,1,2,3 or 4 (see 
right hand side of figure 11). Moreover, as 

n 

dBW T (G) 1 \ (J BW r (G)i = {A g C : |A| = 4} 



i=2 



then theorem 



5.2 



implies that 7Zy\r Vl y(N) has a strictly smaller Gershgorin-type 
region than does N which can be seen in figure 12. Considering the fact that n 
may be quite large this example is intended to illustrate that eigenvalues estimates 
can be improved with little effort if some simple structural feature (s) of the graph 
are known. 



7. Concluding Remarks 

A considerable amount of work has gone into the study of matrices with complex 
valued entries (e.g. [T5J|T71I1])- It is quite possible that there are many more results 
in this area, besides those contained in section 3, which can be extended to the class 
of matrices W nx ™. However, such results are not the focus of this paper. 

The main results of this paper demonstrate that isospectral graph reductions 
can be used to improve each of the classical eigenvalue estimates of Gershgorin, 
Brauer, and Brualdi. Moreover, these graph reductions are general enough that 
this process can be applied to almost any graph G G G" and in particular to 
any graph with complex valued weights. Hence, the eigenvalue estimates for most 
matrices in W" xn and all matrices in C nxn can be improved via our isospectral 
reduction process. Additionally, this process is sufficiently flexible to improve such 
eigenvalue estimates to whatever degree is desired. 

Furthermore, graph reductions computationally do not require much effort es- 
pecially if some particular structural feature of the graph is known. In fact, the 
number of calculations required by such estimates may even be reduced by our pro- 
cedure since nontrivial reductions typically produce fewer regions used to estimate 
the graph's spectrum. These properties indicate the potential usefullness of graph 
reductions with respect to applications and particularly to dynamical networks. 
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Moreover, this paper also raises new questions related to graph reductions and 
eigenvalue estimates. For instance, what algorithms related to choosing structural 
sets and sequences of structural sets can be developed to improve the speed or 
accuracy of such estimates. 
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